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Abstract 

We extend the threshold resummation exponents G N in Mellin-A^ space to the fourth logarithmic 
(N 3 LL) order collecting the terms a^(a s lnA^)" to all orders in the strong coupling constant oc s . 
Comparing the results to our previous three-loop calculations for deep-inelastic scattering (DIS), 
we derive the universal coefficients B q and B g governing the final-state jet functions to order oc s 3 , 
extending the previous quark and gluon results by one and two orders. A curious relation is found 
at second order between these quantities, the splitting functions and the large-angle soft emissions 
in Drell-Yan type processes. We study the numerical effect of the N 3 LL corrections using both the 
fully exponentiated form and the expansion of the coefficient function in towers of logarithms. 



1 Introduction 



Coefficient functions, or partonic cross sections, form the backbone of perturbative QCD. These 
quantities are defined in terms of power expansions in the strong coupling constant oc s . In general, 
only a few terms in this expansion can be calculated. It is however possible, and necessary, to 
resum the dominant contributions to all orders in a s close to exceptional kinematic points. Close 
to threshold, for example, where real emissions are kinematically suppressed, the resummation 
takes the form of an exponentiation in Mellin-/V space [1-4], with the moments N defined with 
respect to the appropriate scaling variable, like Bjorken-x in deep-inelastic scattering (DIS) and 
x T = 2p T l \J~S in direct photon and inclusive hadron production. 

The resummation exponents are given by integrals over functions in turn defined by a power 
series in oc s . Besides by dedicated calculations, the corresponding expansion coefficients can be 
obtained by expanding the exponentials and comparing to the results of fixed-order calculations. 
Hence progress in the latter sector also facilitates improved resummation predictions. At present 
the next-to-leading order (NLO) is the standard approximation for many important observables, 
facilitating a resummation with next-to-leading logarithmic (NLL) accuracy. For recent introduc- 
tory overviews see, for instance, Refs. [5,6]. The next-to-next-to-leading order (NNLO) correc- 
tions have been completed so far only for the coefficient functions for inclusive lepton-proton 
DIS [7-11], the Drell-Yan process [12-14] and the related Higgs boson production [13, 15-17] 
in proton-proton collisions. Consequently, the threshold resummation has been carried out at the 
next-to-next- to-leading logarithmic (NNLL) accuracy only for these processes [18, 19]. 

Recently we have computed the complete three-loop coefficient functions for inclusive photon- 
exchange DIS [1 1, 20]. Moreover, in the course of the calculation of the third-order splitting func- 
tions governing the NNLO evolution of the parton distributions [21, 22], we have also computed 
DIS by exchange of a scalar directly coupling only to gluons. Together these results enable us 
to extend two more universal functions entering the resummation exponents, the quark and gluon 
jet functions B q and B g collecting final-state collinear emissions, to the third order in cc s . In fact, 
already the second-order coefficient for B g represents a new result, relevant for future NNLL re- 
summations of processes with final-state gluons at the Born level. Making use also of the results 
of Refs. [23, 24] we can furthermore effectively, i.e., up to the small contribution of the four-loop 
cusp anomalous dimension, extend the threshold resummation for inclusive DIS to an unprece- 
dented next-to-next-to-next-to-leading logarithmic (N 3 LL) accuracy. 

The remainder of this article is organized as follows: after recalling the general structure of 
the resummation exponents in Section 2, we extend the required integrations in Section 3 to the 
fourth logarithmic (N 3 LL) order. In Section 4 we determine the relevant expansion coefficients 
by comparison to our three-loop results for DIS and illustrate the numerical effect of the N 3 LL 
contributions to the resummation exponent. In Section 5 we present the resulting predictions for the 
leading seven large-* terms of the four-loop coefficient function and discuss higher-order effects 
in terms of the expansion in towers of threshold logarithms. Our results are briefly summarized in 
Section 6. Some basic relations for the integrations of Section 2 can be found in the Appendix. 
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2 The general structure 



For processes with only one colour structure at the Born level, the resummed Mellin-space coeffi- 
cient functions C N (defined in the MS scheme) are given by a single exponential [1,2] 

C N (Q 2 )/C» (Q 2 ) = g (Q 2 )-oMG N (Q 2 )} + O (ATi hi" AO . (2.1) 

Here C^ Q denotes the lowest-order coefficient function for the process under consideration, e.g., 
C^q = 1 for DIS. The prefactor g Q collects, order by order in the strong coupling constant a s , 
all AMndependent contributions. The exponent G N contains terms of the form ln^A^ to all orders 
in oc s . Besides the physical hard scale Q 2 (= — q 2 in DIS, with q the four-momentum of the 
exchanged gauge boson), both functions also depend on the renormalization scale fj r and the mass- 
factorization scale /j jr. The reference to these scales will be often suppressed for brevity. 

The exponential in Eq. (12.11) is build up from universal radiative factors for each initial- and 
final-state parton p, A p and 7 p , together with a process -dependent contribution A mt . For exam- 
ple, the resummation exponents for inclusive deep-inelastic scattering, Drell-Yan (DY) lepton-pair 
production and direct photon production via qq — > gy and qg — > qy [25] take the form 

G£ IS = lnA q + ln/ q + lnAgl S) 
g£ y = 2 InAq + hiApy , 

GJLcy = lnA a + lnA b + ln/ c + lnA^ ■ (2.2) 

The radiation factors are given by integrals over functions of the running coupling. Specifically, 
the effects of collinear soft-gluon radiation off an initial-state parton p = q,g are collected by 

rl t^-I _ 1 riy-zfQ 1 da 2 
lnA p (2 2 ,^ ) = / d z / _«_ Ap(0ii(9 2 )) . (23) 

Jo i—Z Jfjj q 
Collinear emissions from an 'unobserved' final-state parton lead to the so-called jet function, 



ln7 p (2 2 )= / dz-. / ^A p (a s (q 2 ))+B p (a s ([l-z]Q 2 )) 

J0 t Z J (l—z)Q Q 



(2.4) 



Finally the process-dependent contributions from large-angle soft gluons are resummed by 

lnA int (£ 2 ) = f l dz z ^--Z±D{as{[\ ~z] 2 Q 2 )) . (2.5) 
Jo l-z 

The functions g Q in Eq. (12.11) and A p , B p and D in Eqs. (|2.3I) - (|2.5I) are given by the expansions 

= L F i±^ L^s, (2-6) 

l=k ™ l=k 

where Iq = with g 00 = 1 for F = g , and Iq = 1 else. Here we have also taken the opportunity 
to specify the reduced coupling a s employed for the rest of this article. The extent to which these 
functions are known defines the accuracy to which the threshold logarithms can be resummed. 
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The situation for inclusive DIS is actually simpler than indicated in Eq. (|2.2I) . as the function 
(12.51) is found to vanish to all orders in a s [23,24], 

Df 18 = , AgJ s = 1 . (2.7) 

On the other hand, the resummation of processes with four or more partons at the Born level, 
like inclusive hadron production in pp collisions [26], is more complicated than Eq. (12.11 ) due to 
colour interferences and correlations in the large-angle soft gluon emissions [3,4]. However, the 
process-independent functions A p and B p retain their relevance also for such cases. 



3 The resummation exponent to fourth logarithmic order 



After performing the integrations in Eqs. (12.31) - (12.51) . the function G N in Eq. (12.21) takes the form 

G N (Q 2 ) = \nN- gl (l)+g 2 (k)+a s g 3 (X)+a 2 g 4 (X) + ... , (3.1) 

where X = $oa s \nN. For the actual computation of the functions g- it is convenient to employ the 
following representation for the scale dependence of a s up to N 3 LO : 



a s (q 2 ) 



+ 



4|i lnX +4 

X 2 Po X 3 



P 2 



i(ln 2 X 



hiX-l+X) + ^(l-X) 
Po 



+ ^(1 



5 n a 1 1 - 

2(1 -X)\nX + -\n 2 X-\n 3 X--+X--X' 
v 1 2 2 2 



X 2 ) + ^^(2X\nX -3lnX -X(l -X)) 



P 2 , 



+ 0(a s 5 ) (3.2) 



with a s = a s (/u 2 ) and the abbreviation X = 1 + a s $oln(q 2 / /u 2 ) . The terms up to the n-th order in 
a s in Eq. (13.21) contribute to g n in Eq. (13.11) . Thus the calculation of g 4 requires the highest known 
coefficient of the beta function of QCD, p*3 [27, 28]. 

Generalizing the approach of Ref. [18], the functions g ; -(X) can be obtained using well-known 
methods for Mellin transforms based on properties of harmonic sums and harmonic polyloga- 
rithms [29,30] in addition to algorithms for the evaluation of nested sums [31]. The basic relations 
for this approach, suitable for the evaluation of Eq. (13.11 ) to any accuracy, are presented in the 
Appendix. As a check we have also carried out the integrations along the lines of Ref. [19]. 

For the convenience of the reader, we first recall the known results for g^ g 2 and g 3 [1, 2, 
18, 19]. For brevity suppressing factors of Po (see below) and using the short-hand notations 
L qr = \n(Q 2 //j 2 ) and Lf r = \n(/j 2 /iu 2 ) , these functions can be written as 

Ai(l-ln(l-X)+^ 1 ln(l-X)) , 



.DIS/ 



S2 DIS W 



(A 1 p 1 -A 2 )(^ + ln(l-^)) + iA 1 p 1 ln 2 (l-X) 
- (Ai y e - B X ) ln( 1 - X) + L qr A 1 ln(l - X) + L fr Ai 7, 



(3.3) 



(3.4) 
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S 3 DIS W = ^(A 1 p 2 -A 1 (3f+A 2 (3 1 -A 3 )(l+X- T 4 



+ (AiPiy e +A 2 |3i-5i|3i) 1 



1 ln(l-X) 



- (AxPj + ^Aj (y 2 + + A 2Ye - B lJe -B 2 )(l- j^) 
+ L qt [(Aij e -A1P1 +A 2 -Bi) (l - y-^) +A1P1 (^y£^ 
+ L t A 2 l - L\}-A X (l - - L| ^AiX . 



(3.5) 



The dependence on Po is recovered by A^ — > A^/Pq, ^ - ^ #fc/Po% P^ > P^/Po +1 an< ^ rnultiplication 
of g3 by Po- In the same notation the new function #4 (to be multiplied by Pq) is given by 



1 3 ln 3 (l-X) 1 2 2 2 ln 2 (l-X) 1 

" AlPl +2 (AlPlYe+A2,3l " jBl|3l) (1-A,)2 + 2 



'A 1 p 1 (y e 2 + ^ 2 )+A 2 pf-2A 2 p 1 Y e -A 3 pi + 25^ + 25^ 



+ ^(A 1 p3-A 1 p 1 p 2 
In(l-X) 



- (AjP? -AjPjfe)^^ + (±AiP? - AiPiP 2 + ^3) ln(l - 1) + (Atf? 
-A 1 p 1 p 2 -A 1 p 2 y e +A 1 p 2 y e -A 2 p 2 +A 2 p 2 +5 1 p 2 -5 1 p 2 )(i- T l I 

5 1 5 1 

+A 2 piy e - A 2 (y 2 + £ 2 ) - z A 2 pf + -A 2 p 2 + -A 3 pi - A 3 y e - -A 4 -B2P1 

o J o i 

1 



+ 5 1 (y 2 + ^ 2 )+25 2 y e +5 3 )(l-^-i x ^) + i(A 1 p3-2A 1 p 1 p 2 +A 1 p 3 
+A 2 p 2 -A 2 p 2 +A 3 p! _A 4 )X + L qr [(Ajp 2 -A1P2) " T3X + ^ (1^)2 
+ (^Ai (y 2 + £ 2 ) - ^A 2 p! +A 2 y e + l -A, -B lJe -B 2 ) (l 



(!-*,)< 



ln(l-X) 1 R2 ln 2 (l-X) 
ln(l-A,) 



+ (AiPiYe +A 2 pj -^lPO ri^z - ^AiP 2 



+ «AiP r 



+A 2 -B 1 )(] 
+ L ir A^-Ll(A 2 + l -A l $ i y + Li l -A l 'k. 



(l-X) 2 J 2 IKI (1-X) 2 



1 



2 (AlYe 
1 



+ L V6 Al l 1 (l-^) 2 



(3.6) 



The results g ; DY for the Drell-Yan process and, with slightly different coefficients A, and £>,, Higgs 
production via gluon-gluon fusion, are related to Eqs. (13.31) - (|3.6I) as follows: the function cor- 
responding to Eq. (|3.3I) reads gP Y (?i) = 2gP IS (2?i), while the functions g 2 Y ,gf Y and are 



obtained from Eqs. (13.41) - (13.61) by replacing X — ► 2X everywhere and substituting Bi — > Di/2 in 
all terms. Finally, the constants have to be changed according to y e — > 2y e and — > 2 n £ n . The 
generalization of Eqs. (13.31) - (13.61) to other processes involving Eqs. (12.31) - (12.51) is obvious. 

The functions gi(X) collecting the leading logarithms L(a s L) k depend on A\ only and are finite 
for all X. The N" _1 LL contributions g n >i(X) to Eq. (|3.1I) including A n , 5„ i and D n -\, on the other 
hand, exhibit Landau poles at the moments N = exp[l / (pV*s)] f° r DIS and N = exp[l / (2pV*,y)] 
for the DY case (and at both values in general). As obvious from Eq. (13.31) - (13.61) the strength of 
these singularities increases with the logarithmic order, reaching [1 — (2)X]~ 2 at the N 3 LL level. 



4 Resummation coefficients and numerical stability 

Since observables are independent, order by order in a s , of the factorization scale, the functions 
A p in Eqs. (12.31) and (12.41) are given by the large-Af coefficients of the diagonal splitting functions 
for n r = /J f , 

Pp P (a s ) = -A p (a s )ln/V+o(l) , (4.1) 

which in turn are identical to the anomalous dimension of a Wilson line with a cusp [32]. The first 
and second order coefficients have been known for a long time, the third order has been recently 
completed by us. The expansion coefficients (12.61) for the quark case read [21, 33] 



Aq,i = 4C> 

Aq,2 = 8C/r 



l q,3 



160 



6 l-^ 2 \C A --n f 



2 ,245 67 11 11 2 
L * I — 7T + -r 4)3 + T 4>2 



24 



+ C F n f 



h2C3 

24 hi 



'f 



209 10 „ 7 , \ 2 ( 1 



(4.2) 



Here denotes the number of effectively massless quark flavours, and C> and Ca are the usual 
colour factors, with C> = 4/3 and Ca = 3 in QCD. The gluonic quantities are given by 

Ag (i = Ca/CfA^ . (4.3) 

The perturbative expansion of A p is very benign. For ny = 4, for example, Eqs. (14.21) lead to 

A q (a s ) = 0.4244 a s (1 + 0.6381 a s + 0.5 100 a 2 + ...). (4.4) 

Consequently, already the effect of A3 on the resummed coefficient functions is very small [18, 19], 
and a simple estimate suffices for the presently unknown fourth-order coefficients A4 entering g 4 . 
With the [0/2] results differing by less than 10%, we will employ the [1/1] Pade approximants 



A q , 4 w 7849, 4313, 1553 for n f = 3, 4, 5 , 
corresponding to an estimate of + 0.4075 for the next term in Eq. (14.41) . 



(4.5) 
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For the determination of the coefficients 5, we also need the constant-TV piece, g in Eq. (12.11) . 
for inclusive DIS. The corresponding expansion coefficients g Qk can be obtained by Mellin invert- 
ing the +-distribution and 8(1 — x) parts of the k-th order coefficient function c 2 „. Again using the 
expansion parameter a s = a s / (47l), the presently known terms [7, 10, 20] are given by 

goF = C F (-9-2^ + 2y 2 + 3y e ) , (4.6) 

sg 8 = c2(^-^y e -^ + 6y e 3 + 2y e 4 + ^ 2 -18y^-4fe 

rrr nA y ^ ^ ( 5465 3155 367 2 

-66^3 + 24y^ 3 + f£\ + C A C F [—^ + ~^le + j^J 2 

22 3 1139. 22 r A 2r 464. An „ 51 . 2 
+y t - - -jle^i - 4y 2 ^ + ^3 ~ 40y e k + y $ 

„ {451 247 29 2 4 3 85 r 4 r A r \ 
+ Cpn f {~36 ~ "27 Ye ~ ~9 le ~ 9 le + ~9 C2 + 3 Y ^ 2 + 9^ 3 J (4 " 7) 



and 



-Dis C 3f 7255 1001 187 2 93 3 , ^ 4 , 6197 

579 o i/i o 

+ — y^2 + 66y 2 ^ - 36y^ 2 - 4y^ 2 - 41 lfe - 346y^ 3 - 60y 2 C 3 

1 791 R 

+ 48y^ 3 - — Q + 84y^ 2 2 + -y 2 ^ 2 + 556^ 3 - 80y e ^ 3 + 1384C 5 

^ r 8144 r3 176 ^2^ ^^2/^9161 16981 5563 2 
-240y^ 5 + —g " + QC 2 ^_ - — y, - — y 2 

8425 3 433 4 44 , 191545,, 28495 . 592 2y n 4 „ 
+ -j^Je + ~y% + y Y e + 1Q8 £2 - ^-T^2 - — Y,%2 - 8y^ 2 

284 ,„ 49346. . 640 2 . 0rt ,„ 11419. 2 

- -g-Jeti - +752y e ^ 3 + — y 2 ^ 3 - 80y^ 3 + 

23098 536 „,\ 2 / 1909753 599375 50689 2 

-C2 3 + -^ 3 2 + c 2 A c F — — — + -^-y e + —^y 2 



315 ^ 3 ^ A r V 1944 729 ' 162 



, 2 / 341 2003 83 2 683 3 70 4 8 5 10733, 



+ c h {-^e+ W Y ' + igY. - 2ftf - - 9 Y * " 

2177 „ 112 ly 32 ,„ 10766,, 20 „ 8 2 „ % y2 
+ -zj-V'b + + Y 1 ^ 2 + ~2J~^ 3 ~ "9~ Y ^ 3 + 9 Ye ^ 3 ~ 3 Ye ^ 2 
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10802 „ 2 40 



784 



135 S-f fcC3-^C« 



^ ^ / 142883 
+ C A C F n / I - 



160906 



7531 2 
Ye 



1552 3 _ 44 4 33331 

81 ,e ~JT le ~T7 le + ~8]~^ 2 + 
21418.. 1976 „ lr 164 y2 



486 
5264 



81 

+ C F Hf 

116 



27 

9517 8714 



135 



729 

64 
~9 



"Ye 



56 7„ 16 a„ 

Y^2 + ^rY e %2 + ^Y^2 



81 

128 „ 2 64 „ „ 8„ 



+ ■ 



470 



486 729 



"Ye + ^rrY e + 



116 



, 80^ 64 „ 
27 Y^2 + ^3 + ^3 



81 ,e 

292 



81~ Ye + 27 Ye 



2110 



81 



135 



+ =-^//n (64 + 160^2 
n r 



224 



32, 
T 



^3 - 



1280 



(4.8) 



Note the new flavour structure /7n [20,34] in g 03 . This contribution, introducing the colour factor 
d abc dabc / n c , for the first time leads to a difference between the flavour-singlet and non-singlet 
coefficient functions for the photon-exchange structure function F 2 in the soft-gluon limit, with 
flff = 3(e) and flf^ = (e) 2 / (e 2 ) , where (e ) represents the average of the charge e k for the active 
quark flavours, (e k ) = nj L;£i e \ ■ Correspondingly, the large-Af coefficient functions for Z- and 
W -exchange DIS will differ from each other and from F 2 ' m ' at order aj. Based on the size of the 
flu term in Eq. (14.81) . however, we expect these differences to be numerically insignificant. 

Now the coefficients B„ ^ entering the jet function (12.41 ) can be derived successively from the 

In N terms of the /c-loop DIS coefficient functions c 2 . Expansion of Eqs. (13.11 ) and (13.31 ) - (13.61 ) 
in powers of a s yields 



'2,q 

,(2) 
2,q 

,(3) 
2,q 



InN 



InN 



InN 



AlJe-Bi 

^Ajpo (Je + C2) +A 2 Je-Bi%j e -B 2 +g i(Aij e -Bi) 
^iPo(Ye + 3Y^2 + 2^ 3 ) + \a$i (J 2 + + A 2 Ule + & 

'Ye) 



+ A 3 y e -Bi (p lYe + (3^2 + %T e ) ~ 25 2 (3 Y e - B 3 
+ g 2(Aa e -Bi)+g 0l (^Ai ^ (j 2 + i; 2 ) + A 2 j e -Bi%y e -B 2 S ) . (4.9) 

/ (k) 

The coefficients of In N , 2 < I < 2k in c 2 L on the other hand, are completely fixed by lower-order 
resummation coefficients, thus providing an explicit /c-loop check of the exponentiation formula. 
Comparison of the relations (14.91) with the corresponding results from the fixed-order calculations 
of Refs. [7, 10,20], using Eqs. fQ . (1431) and d%7> . leads to 



#q,l 
#q,2 



-3C 



+ C F n 



-- + 12^ 2 -24C 3 



247 8 

~in~ 3 



+ C F C A 



3155 44 r ' 
^ + T ^2 + 40C 3 



(4.10) 



(4.11) 
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B 



q.3 



29 

— -I8C2-68C3 



288 „ 2 



C 2 z + 32^3 + 240^ 



+ C A C F 



+ C A C F 



712 

-46 + 287^2-^-^3 



272 



^1-16^2^3-120^5 



599375 32126 „ 21032 „ 
+ ~ C2 + Cs 



729 



81 



+ Cpn 



5501 



50^2 + f 



+ C A C F n f 



54 

160906 9920 



729 



81 



27 
+ C F nj 
776 



652 



176 



^3-232^ 



^3 + 



8714 232 , 32 „ 

1 C2 ^3 

729 27 b 27 bJ 



208 2 

4>2 



15 



(4.12) 



Eq. (14.101) is, of course, a well-known result [1,2]. Eq. (14.1 II) has been derived by us before [35], 
establishing D® 1S = from the n f \n 2 N term at three loops. For our new result (14.121) , on the other 
hand, we have to rely on the subsequent all-order proofs of Eq. (12.71 ) in Refs. [23, 24]. The QCD 
expansion of B q analogous to Eq. (14.41) appears far less stable than that for A„, 



5 q (a s ,n, = 4) = -0.3183 a s (1 - 1.227 a s - 3.405a 2 + ...) 



(4.13) 



The ingredients for the resummation of inclusive DIS are now complete, and in the left part of 
Fig.[T]we show the corresponding LL, NLL, N 2 LL and N 3 LL approximations to the exponent (13.11 ) 
resulting, for a s = 0.2 and three flavours, from Eqs. (T33T) - (TJB . (1421) . (1431) and (14.101) - (14.121) . 
For these parameters the expansion (13.11 ) is stable in the /V-range shown in the figure. For example, 
the relative N 3 LL corrections amount to 2% at N = 10 (k = 0.33) and 4% at N = 40 (k = 0.53), 
whereas the corresponding N 2 LL figures read 9% and 12%. The large third-order contribution to 
B q actually stabilizes g 4 (A,): for B q j = the N 3 LL term at /V = 40 would instead reach 12%, i.e., 
the size of the previous order. The effect of both A q 4 and P3, on the other hand, is very small, as 
their respective nullification would change the result even at N = 40 by only 0.6% and 0.1%. 

In the right part of Fig. 1 and in Fig. 2 the exponentiated results are convoluted with the typical 
input shape xf = x 5 (l — x) 3 for a couple of values for a s and n^. The Mellin inversion is in 
principle ambiguous due to the Landau poles briefly addressed at the end of Section 3. We employ 
the standard 'minimal prescription' (thus adopting the usual fixed-order contour) of Ref. [36], to 
which the reader is referred for a detailed discussion. For total soft-gluon enhancements up to 
almost an order of magnitude, as shown in the figures, the resulting N 3 LL corrections remain far 
smaller than their N 2 LL counterparts and amount to less than 10% even for oc s = 0.3. Note that 
the dependence on is larger than the effect of g 4 . Thus, at this level of accuracy, a reliable 
understanding of heavy-quark mass effects is called for also in the limit x — > 1. 

The gluonic coefficients corresponding to Eqs. (14.101) - (14.121) can be obtained in the same 
manner from DIS by exchange of a scalar (|) with a pointlike coupling to gluons, like the Higgs 

(k) 

boson in limit of a heavy top quark. We have derived the corresponding coefficient functions cl ' 
up to k = 3 already during the calculations for Ref. [22], as a process of this type is required to 
access the lower row of the flavour-singlet splitting function matrix. Comparing those results to 
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Figure 1 : Left: the LL, NLL, N 2 LL and N 3 LL approximations for the resummation exponent for 
standard DIS. Right: the convolutions of the exponentiated results with a typical input shape. 
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Figure 2: As the right part of the previous figure, but for a different value of n* (left) and a s (right). 



Eqs. for cf] yields 



11 ^ 2 



Po, 



B 



g,3 



Ci 



611 88 r , ' 
■-^- + ^^2 + 16^3 



1492081 60875 
21014 



1458 
498329 



+ Car 



1458 
' 8579 
_^4~ 
48829 



■ + 



81 

16^2 
716 



^2 + 



13796 



428 



27 



T ^2 



20 2 
+ 2C F n f -—n f , 



(4.14) 
(4.15) 



27 
296 



2596 



128 



C2C3- 112 ?s 



r a. 568 r 2 

<=3 + TF" 4,2 



15 



CpHf 



832 



r 32 2 



+ C>R 



47 32 ' 



176 



+ m 



200 
243 



8 r 
9^2 



(4.16) 



1458 27 " 27 

where Eqs. (14.151) and (14.161) are new results. For = 4 the numerical expansion of B g reads 

fl g (oc s ) = -0.6631 oc s (1 - 0.7651 cc s - 2.696 c£ + ...) , (4.17) 

exhibiting an enhanced third order correction similar to that of B q in Eq. (14. 13b . 

The gluonic threshold resummation resulting from Eqs. (|4.3I) and (14.141) - (14.161) is illustrated 
in Fig. 3 using the practically irrelevant scalar-exchange process, with same parameters as in Fig. 1 
for direct comparison. The soft and collinear radiation effects are much larger here due to the larger 
colour charge of the gluons, but the qualitative pattern is rather similar to 'normal' inclusive DIS. 

As mentioned above, the (closely related) Drell-Yan process and Higgs boson production via 
gluon-gluon fusion presently represent the only other processes for which the NNLL threshold 
resummation is known, with D\ = and [18, 19] 
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Extending a result given in Ref. [18], we notice the following conspicuous relation between these 
coefficients and 5 p 2 : 
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(4.19) 



where denotes the coefficients of 8(1 —x) in the diagonal two-loop splitting functions, and 

the colour structures on the right-hand sides are those of A p j and B v j = — P g , multiplied by Po- 
Note especially the non-trivial cancellation of all ^-function terms between the three contributions 
on the left-hand sides of Eqs. (14.191 ) and the vanishing of the right-hand sides for Po — > 0. 
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5 Fourth-order predictions and tower expansion 

Another manner to organize the all-order information encoded in Eqs. (|2.1I) - (12.51) is to re-expand 
the exponential, 

2k 

C N (Q 2 )/C? (Q 2 ) = ! + £«* ZcuhP-'+'N, (5.1) 

k=l 1=1 

and retain only those terms in the second sum which are completely fixed by the available infor- 
mation on the expansion coefficients in Eq. (|2.6I) . Using the notation 

oo 

gi$) = Y.Sik^ (5.2) 

k=l 

for the expansion of Eqs. (13.31) - (13.61) together with Eq. (12.61) for g Q , the quantities cu in Eq. (15.11) 
receive contributions from the following coefficients: 
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The complete relations for the first four terms c kl ... c k4 can be found in Ref. [37] in a slightly 
different notation, g 3l — > g 32 (i.e., the second index denoting the total power of a s in Eq. (|3.1I) ). 
Note that the quantities c kl vanish factorially for k — ► °° and fixed /. 

Taking into account Eq. (12.71) and considering the coefficient A,- as either known or irrelevant, 
the function g, for inclusive DIS is completely specified by its leading term, obtained by matching 
to the z'-th order calculation of the coefficient functions as in Eqs. (14.91 ) . The same holds for other 
processes, at least for cases like Eqs. (12.21 ), once the required coefficients 5 ; are known from DIS. 

The leading three towers of logarithms, c kl for any k and / = 1 , 2, 3 , are fixed by a one-loop 
calculation (providing g a for i = 0, 1,2) together with the NLL resummation (adding g lk and g 2k 
for k > 2). This is the status for many important observables [5, 6]. Correspondingly, a two-loop 
computation of the process under consideration specifies g 31 and g 02 and hence fixes, together with 
the NNLL resummation, also the next two towers. This is the accuracy reached for the Drell-Yan 
process and Higgs production [18, 19]. Finally a three-loop computation combined with the N 3 LL 
resummation fixes the first seven towers, c kl for I = 1 , . . . , 7 . With the results of Ref. [20] and 
Sections 3 and 4, we have now reached this point for the structure function F2 in DIS. 



The resulting four-loop predictions, in x-space expressed in terms of the coefficients of the 
+-distributions <D\ = [(1 — x)~ l ln(l — *)]+, for the six highest terms read 
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The seventh term with ©1 is not exactly known, since the fourth-order contribution to A q has not 
been computed so far. Inserting the numerical values for the ^-functions and the QCD colour 
factors, including d abc d a b c /n c = 5n f /18, the resummation prediction is given by 
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-286702 + 642 19.0 /if - 2019.24 nj + 2.0166 nj - 63.402 fl n n f +A qA . (5.10) 



As mentioned above, numerically insignificant are both the uncertainty due to A q 4 (estimated in 
Eq. (14.51) ) and the singlet / non-singlet difference introduced by the //] 1 contribution of Eq. (14.81) . It 
is also interesting to note that the fourth coefficient P3 of the beta function [27,28] with its quartic 
group invariants d F bcd and d A only enters the eighth tower, starting at the fifth order in a s . 



13 





K 1 




KJ 






/CO 


K 1 


1 
1 


Z.OOOO / 


/.U /o4 












Z 


.5. J j J jo 


ZD.ZOJ4 


4U. /oU 


— O/.Ij 








1 
J 




A A OO 1 f\ 
44.yzlU 


zJo.ooj 


4 /U.oZ 


— DZU.J 






A 

4 


Z. 1U /UU 


/I C 

4o. /uyu 


4/ /.oj4 


O/IOO /I A 
z4zy.4u 


CO/IA A 

JZ4U.U 


1 CO/I 

— 18Z4 


OAO 1 


c 

J 


1 1T3T3 


JO.-5ZJ4 


CO 1 C1Q 


jUIj.Io 


zj 1 jU.j 


A QA CO 

4o4oZ 


1 1 O^Q 


D 


A A QO/I /I 


ZJ.JOl / 


CAC QTO 

jUj.y /z 


D4jz.yj 


jzlzy. / 


ZZjJZU 


D / jUIz 


n 

I 


A 1 QAO/I 

u. lyuzo 


1 1 Q<O0 


l/IA QC/I 

J4U.yj4 


com z:i 


DouUz.y 


/I C<~7 1 O 
4oj / 1Z 


OA QA QA 1 
Z4y4o41 


8 


0.06342 


5.0463 


186.822 


4244.45 


65550.0 


668223 


4979993 


9 


0.01879 


1.8583 


86.041 


2467.72 


48805.8 


666670 


6718531 


10 


0.00501 


0.6028 


34.118 


1204.34 


29604.7 


517490 


6747332 



Table 1 : Numerical values of the four- flavour coefficients c M in Eq. (15.11 ) for the quark coefficient 
function in DIS. The first six columns are exact up to the numerical truncation, and the same for 
F\, F% and F3. The seventh column refers to F\ and F2 for the photon-exchange flavour- singlet 
case, flu = 1/10 [27], and uses the estimate (14.51) for the four-loop cusp anomalous dimension. 

The numerical values of the TV-space coefficients c kl in Eq. (15.11) are presented in Table 1 
for / < 7 and k < 10. Recall that also these coefficients refer to an expansion in a s = GC S / (4rc). 
Whatever the normalization of the expansion parameter, however, the coefficients in each column 
(tower) finally vanish for k — > °°, as mentioned below Eq. (15.31) . Thus the series (15.11) converges 
at all /V 7^ for any finite number of towers Z max , i.e., with the upper limit in the second sum 
replaced by / max - The Mellin inversion of the product with the parton distributions / is therefore 
well-defined, in contrast to the fully exponentiated result discussed above. 

Before we turn to the higher-order predictions, it is instructive to compare the approximations 
by the leading large-* and large-TV terms to the completely known two- and three-loop coefficient 
functions [7, 10, 20]. This is done in Fig. 4 for the successive approximations in terms of the 
+-distributions defined above Eq. (I5.4I ). The corresponding results for the expansion in powers 
of In TV are presented in Fig. 5. Obviously both expansions reproduce the exact large-* behaviour 
(up to terms not increasing as x— > 1 for the ratios shown in the figures) at order a" once all enhanced 
terms, with k = 0, . . . , 2n— 1 or In' TV with / = 1, . . . , 2n, have been taken into account. The 
x-space expansion, however, would lead to a gross overestimate if only the terms up to k ~ n+1 
were known. The convergence in the TV-space approach, on the other hand, is much smoother, with 
a good approximation already reached after n terms. 

A shndar pattern is found for the fourth-order coefficient function eg iHustrated in the S a m e 
manner in the left part of Fig. 6: the expansion in decreasing powers of In TV stabilizes after the 
fourth term. Based on these results and the higher-order coefficients shown in Table 1, we expect 
that the first / logarithms should provide a good estimate up to about the l-th order in a s , but 
severely underestimate the effect of the coefficient functions of much higher orders. Consequently, 
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Figure 4: Convolution of the two-loop (left) and three-loop (right) contributions to the DIS coeffi- 
cient functions C2, q with a typical input shape. Shown are the full results [7, 10,20] and the large-.* 
expansion by successively including the +-distributions T>k, respectively starting with Dj and ©5. 
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Figure 5: As the previous figure, but for large-TY expansion in terms of decreasing powers of In TV . 
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the tower expansion should underestimate the corrections towards x—*l, where more and more 
orders become relevant. This is exactly the pattern shown in the right part of Fig. 6, where the 
predictions of all effects beyond order oc s 3 are compared between the tower expansion and the full 
exponentiation (for the latter again using a 'minimal-prescription' contour [36]). Both approaches 
agree very well, for the chosen input parameters, at x < 0.93, but start to diverge at x > 0.95 where 
the exponentiation is also intrinsically more stable. 




0.8 0.85 0.9 0.95 0.8 0.85 0.9 0.95 

X X 



Figure 6: Left: the successive approximations of the four-loop coefficient function c 2 q by the 
large-Af terms specified in Table 1, illustrated by the convolution with a typical quark distribution. 
Right: corresponding results for the effect of all orders beyond a 3 as obtained from the tower 
expansion with up to seven towers and from the exponentiation up to N 3 LL accuracy. 



6 Summary 

We have extended the threshold resummation exponents [1,2, 18, 19] for few-parton processes to 
the fourth logarithmic (N 3 LL) order collecting the terms a^(a s lnN) k to all orders in a s . For our 
reference process, inclusive deep-inelastic scattering (DIS), the N 3 LL contributions are specified 
by two universal expansion parameters: the four-loop cusp anomalous dimension A q 4 and the 
third-order quantity 5 q ,3 which defines the jet function resumming collinear radiation off an un- 
observed final-state quark. The former coefficient has not been computed so far, but can be safely 
expected to have a very small effect of less than 1%. In fact, the perturbative expansion up to 
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A3 [21, 22] does not exhibit enhanced higher-order corrections, and A4 can be estimated by Pade 
approximations. We have calculated the more important second coefficient 5 q 3 by comparison of 
the expanded resummation result to our recent third-order calculation of electromagnetic DIS [20]. 

The perturbative expansion of B q seems to indicate, as far as this can be judged from the first 
three terms, the onset of a factorial enhancement of the higher-order coefficients. However, the 
rather large size of the coefficient B q ?, actually stabilizes the logarithmic expansion of the coeffi- 
cient functions. In fact, the N 3 LL corrections are very small at large scales Q 2 , and even facilitate 
a reliable prediction of the soft-gluon effects at scales as low as Q 2 m 4 GeV 2 (corresponding 
to a s ~ 0.3) down to very small invariant masses W of the hadronic final state in ep — > eX, 
W 2 — m 2 m 0.5 GeV 2 . Thus we expect our results to be useful also for low-scale data analyses 
using parton-hadron duality concepts. 

The threshold resummation can also be employed to predict, order by order in oc s , the leading 
\\\N contributions to the higher-order coefficient functions. At the level of accuracy reached in 
the present article for inclusive DIS, the exponentiation fixes the seven highest terms, with 
21 — 6 < n < 2k, at all orders k > 4 of a s . Already the highest k powers of InN provide a good 
estimate of the soft-gluon enhancement of the fc-loop coefficient functions at least for k < 7, in 
contrast to the (expected, see Ref. [36]) worse behaviour of the corresponding expansion in x-space 
+-distributions. Except very close to threshold, where too many orders in oc s become important, 
the summation of the above seven N- space logarithms to all orders yields a good agreement with 
the exponentiated coefficient function. This agreement further confirms the 'minimal prescription' 
[36] used for defining the in principle ambiguous Mellin inversion of the resummation exponential. 

Besides the standard (gauge-boson exchange) process, we have also considered DIS by ex- 
change of a scalar directly coupling to gluons. By comparison of the resummation to our unpub- 
lished three-loop coefficient function for this process we have derived, for the first time, the second 
and third order contributions to the coefficient B g governing the jet function of a final-state gluon. 
The quantity 5 g 2 will be employed to extend the NNLL resummation to more processes, once the 
corresponding NNLO results required to fix the process-dependent large-angle soft contribution 
become available. Finally we would like to draw attention to the curious relation (14.191 ) which 
connects, for both the quark and gluon channels, the second-order splitting functions, jet functions 
and the two-parton (Drell-Yan) large-angle soft emissions in a non-trivial manner. 
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Appendix 



Here we show some key elements of the calculation of the resummation exponents gt presented in 
Section 3. Useful auxiliary relations (for \x \ < 1) are 

_2 = y r(n-e+iV 

(i-*) B - e " h nn-e) ii ' 1 ■ ; 

tfli-x) /aV i (A _ 2) 



(l_jt)»-e \d e J (l-x)«- 

The Mellin transforms of the +-distributions follow from the results for harmonic polylogarithms 
[30] and are given by 

l 

I II 1 ■ I . ) 1 I 

l-z 



dz^-^ln k (l-z) = (-\) k+l k\S h ^ A (N) , (A.3) 



o t+ i 

where S mir .. jmk (N) denotes the harmonic sums [29]. Eqs. (IA.1I ) - (IA.3I) lead to the master formula 
for the derivation of the functions gu 

fd Z — 1 = 

J Z l-z (l+aln(l-z))"- £ ~ 

-t |^(*^ | o ( '" + r 1 ) (-pomf y (a., 

i+i 

with a = (a s Po)/(\+a s fi \nQ 2 /tf). The double sum in Eq. (IA.4I) can be solved to the desired log- 
arithmic accuracy with the algorithms for the summation of nested sums [31] coded in Form [38]. 
The expansion of the Gamma function in powers of £ for positive integers n reads 

T n\T{\+€) = 1 + e5 lW+ e2 ( 5 MW-^W) + £ 3 (5 1 , u (n)-5 1 , 2 (n) 
-52,i(n)+5 3 H)+£ 4 (Si,LLi(n)-5 u ,2W-Si,2,iW 
+ Si, 3 (n)-S2,i,i(n)+S2,2(n)+S3A(n)-S 4 (n)) + 0(e 5 ) . (A.5) 
Finally, with 0/ y = 1 for i > j and 0/ 7 = else, the sums Si,.„,i(iV) are factorized according to 

= (smy+\Ki-mmsim i - 2 +lKi-m-2)s 3 (N)(s l (N)y- 3 

+ l -i(i-l)(i-2)(i-3)(s 4 (N) + 1 -(S 2 m 2 ysmy- 4 + --- (A.6) 

~ a In' A/ + ie ; - 3 i(^ - l)^ 2 ln^ 2 + |0 r - 4 ^(^ ~ 1) ~ 2)^ 3 In'" 3 A> 

+ie B /(« - 1) (/ -2) (i - 3) fa + ^2 2 ) In- 4 A>+ . . . , (A.7) 
where N = Ne^ e and the algebraic properties of harmonic sums have been used [29]. 
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